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ATTITUDE CONTROL AND DRAG COMPENSATION PROPULSION SYSTEM 
FOR THE GRAVITY PROBE-B SPACECRAFT 

I. INTRODUCTION 


The on-board propulsion system for the Gravity Probe-B spacecraft must provide attitude control 
plus thrust to effectively cancel the atmospheric drag encountered at the relatively low orbit dictated by 
current shuttle capability. A propulsion system was desired which would utilize die boiloff of helium 
gas from the experiment package dewar as propellant. This helium boiloff would normally be vented 
non-propulsively and, thus, in a sense be “wasted.” A system was devised which could reconcile the 
thrust requirements and the available helium boiloff to provide 100 percent drag free control of die 
vehicle. 


Derivation of the necessary thermodynamic equations is shown and a computer code for a rela- 
tively rigorous simulation of the two-phase helium dewar is provided along with details of die delivery 
line pressure loss calculation. These mathematical models were used to estimate the dewar pressure and 
temperature as a function of time during the predicted extremes of thruster demand and dewar heat 
absorption. 


II. SYSTEM REQUIREMENTS AND DESIGN FACTORS 


The requirements for the propulsion system were established by the need to provide precise 
attitude control and to provide an additional thrust vector to effectively cancel die vehicle drag incurred 
by the spacecraft The thrusters are paired back-to-back with clusters of two pairs located at the ends 
of the four solar array spars for a total of 16 individual thrusters. Their location relative to the vehicle 
is illustrated in Figure 1. Each of the thrusters consists of a flow control valve and nozzle and is required 
to produce a maximum thrust of 650 dynes with a thrust tolerance of ±2 percent of the open loop 
commanded value. The desired transient response of the thruster to a step input thrust command is that 
of an underdamped second-order system 


e n } T 

f / f COM ■ 1 + / •*» ty/l-t'-n) 

/lTf 5 

where 

a = cos** K) in radians 
t = time in sec 


F = delivered thrust 




Gravity Probe-B Spacecraft attitude control system. 



^com * command thrust 

<*> n ■ undamped natural frequency radians/sec 

f ■ damping factor 


with the values of natural frequency, « n , and damping factor, f, set at 10 Hz (6Z8 RAD/sec) and 0.707 
respectively. 

The thruster flow control valve is required to respond to thrust command changes as small as 
0.1 dynes in the closed loop mode. This very small thrust tolerance is necessary to maintain the vehicle 

within the acceleration limit of 1(F 10 g (this assumes two thrusters are operating together for drag 
makeup). 

The maximum total thrust required at any one time is ret at 1300 dynes. Although this thrust 
happens to be equal to twice the maximum thrust of two nozzles, as many as eight thrusters could be 
operating simultaneously at less than their maximum values. This maximum thrust is ret by the maxi- 
mum atmospheric density expected during a 90 min orbit The required thrust will statistically vary as 
shown in Figure 2. 
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Figure 2. Orbital variation of required total thrust 
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The constraints within which the helium bofloff thruster system must operate are set by the 
experiment package dewar design, the 1 year life requirement, and the solar heat pickup capability of 
the helium delivery lines. The heat leak into the dewar (which contains «n initial 300 kg load of helium 
in a saturated liquid form) is expected to range from • low of 0.1849 to a high of 0.2326 W. A nominal 
dewar temperature of 1.6°K with an allowable tolerance of ±0. 1*K was considered aa t j sfoctoty to pro* 
vide the superfluid helium II conditions required by the experiment package. The tolar heat picked up 
by the bofloff helium as it cools die dewar thermal shields mid traverses the lines leading to die thruster 
modules at die ends of die solar array spars raises its temperature to between 300° and 375*K. Finally, 
the vehicle solar array and power supply system indicate a power/weight penalty of 1.5 kg/W. Them 
requirements are summarized in Table 1. 


TABLE 1. SYSTEM REQUIREMENTS AND DESIGN FACTORS 


Satellite Life (minimum) 

1 year 

Orbit Period 

90 min 

Dewar Initial Helium Load 

300 kg 

Maximum Thrust per Nozzle 

650 dynes 

Thrust Tolerance Open Loop 

±2 percent of commanded value 

Transient Response Parameters 

Damping Factor* 0.707 
Natural Frequency * 10 Hz 

Sensitivity (minimum thrust increment) Closed Loop 

0. 1 dynes 

Allowable Unbalanced Drag Acceleration 

i<r 10 g 

Maximum Total Thrust Required 

1300 dynes 

Helium Dewar Heat Leak Range 

0.2326 to 0.1849 W 

Nominal Dewar Temperature 

1.6°K 

Allowable Dewar Temperature Range 

1.5 to 1.7*K 

Boiloff Delivered Temperature Range 

300° to 375°K 

Power Penalty (system man per W) 

1.5 kg/W 
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III. HELIUM THRUSTER PERFORMANCE 


The ideal (one-dimensional, frictionless) vacuum sp ec ific impulse of a helium nozzle at an ami 
ratio of four was calculated using the ODE computer program [1 ] for a dumber temperature and pressure 
of 300°K and 3 torr, respectively. This yielded a characteristic velocity of 3566 ft/sec and a vacuum 
specific impulse of 168.9 lbf sec/lbm. However, at file low dumber press ures and small throat diame- 
ters, the throat Reynolds’ number will be quite low (approximately 100) and there will be considerable 
friction loss. Tests conducted at Stanford University (2] on a 20 deg half angle conical nozzle with a 
1.37 mm (0.054 in.) diameter throat, and area ratio and dumber conditions as above, gave vacuum 
specific impulse figures of approximately 130 lbf sec/lbm. This spedfic impulse effidency of approxi- 
mately 0.77 has been substantiated in as-yet unpublished tests at MSFC with a nozzle of similar geometry 
but a larger throat diameter of 5.01 mm (0.196 in.). It therefore appears that a specific impulse of 
130 lbf sec/lbm, or possibly larger, is attainable for the planned propulsion system. However, for the 
system developed in this report, a value of 1 15 lbf sec/lbm was used for file spedfic impulse to provide 
some contingency. Even this conservative value is very good in comparison with alternative propellants, 
particularly when it is available as a waste pa from file dewar boiloff. 


IV. PROPULSION SY8TEM DESCRIPTION 


The main problem to be overcome on the supply side of the propulsion system is in accom- 
modating (with minimum waste) the inequality between the boiloff produced by an arbitrarily varying 
heat leak into the dewar and the on-demand varying flow rate required by file thrusters. Conversion of 
the maximum dewar boiloff due to heat leak (0.2326 W) and the maximum required thrust (1300 dynes 
at 1 15 lbf sec/lbm spedfic impulse) into helium flow rates gives values of 9.75 and 1 1.53 mg/sec, respec- 
tively, and there is no guarantee that the maximum heat leak would occur at the same time as the maxi- 
mum thruster requirement. Obviously, at other times the thruster demand would be less than the avail- 
able boiloff and a surplus of propellant would be available. 

A system which could save or store the excess flow to be used to cover the periods of deficiency 
is most desirable, especially if it can be accomplished without a weight or power penalty. Such a system 
is made possible by exploiting the small, but significant temperature range allowable for the dewar. The 
large spedfic heat capadty, low heat of vaporization, low heat leak, and large helium mass relative to 
the thruster demand flow rate all combine to cause even extremes in the thruster demand and heat leak 
to result in very small dewar temperature changes over long periods of time. This allows variations in 
flow rats and heat leak to average out, thus minimizing any dewar venting or heating. Avenging the 
heat leak extremes gives an average boiloff of 8.75 mg/sec and using Figure 2 gives an orbit average 
demand of approximately 4.44 mg/tec (500 dynes weighted average thrust at 1 15 lbf sec/lbm spedfic 
impulse). Therefore on a per orbit basis there should be ample propellant boiled off. A schematic of 
the system capable of accomplishing this is shown in Figure 3. The dewar temperature would be con- 
trolled within the allowable range of 1.5° to 1.7°K by activation of a heater on the tank at the low 
temperature limit and venting non-propuisively at the high temperature limit The heater would require 
negligible power due to the very small mass flow rates required and the low latent heat of vaporization 
of helium. The non-propulsive vent logic would be combined with the attitude control/drag makeup 
system. When the high temperature limit is approached and venting needed, file controller would select 
a thruster which was not at full thrust and increase its thrust and that of the opposite one of the back- 
to-back pair an equal amount. This would effectively vent the extra helium boiloff without the cost 
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Figure 3. Helium propulsion system schematic. 

and complexity of an additional n on-propulsive vent valve and lines. Although temperature is indicated 
as the control variable, pressure could be substituted as the measured variable since it is related to the 
temperature by the two-phase conditions existing inride the dewar. 

A sketch of the distribution system from the dewar to the thrusters is shown in Figure 4. The 
helium exits the dewar through a porous plug which acts as a phase separator, keeping the liquid helium 
inside the dewar and allowing the vapor to escape. After the helinm leaves the porous plug, it coeds 
the dewar heat shields using the four cooling rings located around the neck tube and exits the dewar 
into a 1.75 in. OD supply line. This line carries it down tile side of the dewar to the support ring (with 
both sides coupled by the cross braces) which is used as a manifold to distribute the helium to the insid- 
of the four solar array span. A cluster of four thrusters is fed from the end of each spar. Solar heating 
of the dewar heat shields, lines and span raises the helium temperature to at least 300°K by the time it 
reaches the thrusters. Use of the support ring and solar array spare to deliver the helium saves the weight 
of additional lines and the large size both reduces pressure loss and provides volume (approximately 1 17 
litere) to reduce pressure transient fluctuations. 
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V. SYSTEM ANALYSIS AND SIMULATION 


To property evaluate the flow averaging effect of allowing the dewar temperature to vary to effec- 
tively store liquid helium along with energy to convert it to vapor on demand, a relatively rigorous 
thermodynamic analysis of two-phase helium storage was completed. This analysis included heat addition 
to the dewar; vapor removed for thruster demand; changes in the temperature, pressure, and density of 
the stored helium; and transition from one phase to the other. Details of the derivation of the equations 
and representation of the thermodynamic properties are covered in Appendix A. 

A pressure loss equation for the helium gas delivery r.ystem is derived in Appendix B and includes 
use of an average density through each individual flow resistance. A constant temperature is assumed to 
exist across each resistance element although the temperature can vary from element to element 
Checking the Knudsen number, Nj^ n , inside the 1 in. OD (0.02 in. wall) dewar tube which is the most 

stringent point in the delivery system; (Le., where the highest Knudsen numbers would b expected); 
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where 

K = Boltzmann constant ■ 1.38 x KT 16 erg/°K 
T * temperature = 213°K 
P = pressure = 3.0 ton- 

5 3 molecular diameter (helium) * 2.18 x IS'*® cm 
r * tube inside radius = 1.22 cm 


verifies die use oi* tne cont’nuum equations since die transition range from continuum to free molecular 
flow is usually ac cepted as 


0.01 <N Rn < 10 


The resistance factors for use hi the delivery system pressure loss calculation are determined in 
Appendix C. Details on line sizes and lengths can also be found in that section along with flow curves 
for the porous plug at die dewar exit. 

A computer simulation based on the equations and factors developed in the three appendices is 
present 4 *! in Appendix D, along with a sample case. The code was used to evaluate die temperature 
changes at the extreme (worst case) conditions and the results are shown in Table 2. The first two cases 
explore the extremes of heat leak and thruster demand for die full dewar with cases three and four the 
same, but with only 5 percent of the helium remaining. Case five was used as the sample case in 
Appendix D to get larger, more obvious changes and has a very small amount of helium and an unreal- 
istically low (i.e.; 0.0) heat leak coupled with maximum demand. All of these cases show diat 150 min 
of extreme conditions can be tolerated before a temperature limit condition is even approached (starting 
from a nominal 1.6°K condition). This means that temperature fluctuations during more than one 90 
min orbit can be averaged out. Pressure loss through the distribution system at die maximum thruster 
demand rate of 0.01 153 g/sec, and at a low dewar pressure of 4.20 torr (1.533°K saturation temperature) 
was found to be 0.385 torr which is considered satisfactory. 


CONCLUSIONS 


A propulsion system utilizing the boiloff helium gas from the experiment package dewar can be 
made feasible by exploiting the small allowable temperature range of the two-phase helium dewar. lids 
small temperature fluctuation allows liquid helium and sufficient heat to vaporize it to be stored in the 
dewar for later ure when a high thruster demand flowrate is required. T* temperature limits of 1.5° 
to 1.7°K allow a sufficient range to average the expected variations in hvat leak and thruster demand over 
more than one 90 min orbit even under wont case conditions. 
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TABLE 2. COMPUTER SIMULATION RESULTS FOR EXTREME CONDITIONS 


Initial Temperature 1.6°K 
Saturation Pressure 5.689 torr 



Conditions 

Remits® 150 min 

Case 

Initial Helium 
Mass 
(grams) 

Demand Flow 
(g/sec) 

V 

Heat Leak 

(W) 

% 

Temperature 

(°K) 

Saturation 

Pressure 

(torr) 

1 

300,000 

0.01153 (Max) 

I 

0.1 849 (Min) 

1.598 

5.650 

2 

300,000 

0.00 (Min) 

0.2326 (Max, 

1.605 

5.800 

3 

15,000 

0.01153 (Max) 

0.1849 (Min) 

1.588 

5.388 

4 

15,000 

0.0 (Min) 

0.2326 (Max) 

1.634 

6.561 

5 

550 

0.01153 (Max) 

0.0 (Below 
Mia) 

1.533 

4.204 
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HELIUM DEWAR THERMODYNAMIC DERIVATION 



First law energy balance (neglecting potential and kinetic energy): 


6Q - 6 W + (u + Py)^ 6 - (u + Py)^ m mt - d U 


or in time rate of change with 8 W and 8 M- n equal to 0. 


qdt-hy ~| dt * d U 

dt I m,* 


J out 


. dU 
, ' h ' ,w= dT 


The dewar internal energy, U, is the sum of the vapor and liquid portions 
U = My Uy + M L u L 

Substituting the following equalities 

M = pV 
u = h - Pv 
P = 1/v 

v v *v D -v L 

hy = h L + hfg 


( 1 ) 
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and rearranging gives an expression for the internal energy in properties dependant on the saturation, T, 
and a time dependent variable, the liquid volume, V^. 


U * V^, (py (hL + hfg) - P) - V L (py (h^ ♦ h^ - P l 


Functionally then, U depends only on dewar temperature and time 


U * fCT.t) 
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in differential form 


dU 


aul au] 
•TJt «J t 


dT 


and the time derivative of the internal energy is 

du = aul aul dr 
dt at J T aTj t dt 


( 2 ) 


Evaluating the time partial using the expression for internal energy above gives 

dul 9V L 

— Jt - - 0>, (h L + hf g ) - h L> — 


An expression for the liquid volume was found as follows (assuming the mass in the tank exceeds py Vp). 
The total mass is the sum of the mass of the vapor and liquid 


Mj = Py Vy + p L V L * py (V D - V L ) + P L V L 


and solving for liquid volume 


V 


L 


Mj - py Vp 
(p L -py) 
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If tiie thruster demand, W, is constant, the existing mats in the dewar at any time is the initial loaded 
mass less the flow rate times the elapsed time 


M T = Mi - W t 

Substitution gives the desired liquid volume relationship 
Mi-Wt-p v V D 

^ 0>l~Pv^ 

and the time partial derivative of the liquid volume is 

av L _ W 
3t {p\^~ Py) 
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giving the time partial derivating for internal energy below 



The temperature partial derivative for the internal energy was found similarly 

aul ap V P aV L 

I = " V L> < h L + bfg) Jf + [^L ' PV> h L * PV h fg] 

ah L ah fe 9P 

+ [py v d + (PL ~ p V) v L ] w + py CV D - V L> - v d ax 

and, from the expression for liquid volume, the temperature partial derivative for liquid volume is 

0V L _ Mj - Wt - p v V D V D dpy 
aT (p l -p v ) 2 (pl-pv) 3T 
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After substitution and simplifying, the temperature partial derivative for internal energy is given below 


au 

3T 




3P 


dhL Py 3hfg &Py ^ Py 


Convi 


3T (PL-Py) 3T 0>L~Py) dT L (PL~Py)J 


Rearranging equation (2) to get the dewar saturation temperature time derivative 

du aul . r p V i . 

dT dr-irJ T ‘-*-< h L^ g )w.-[ 5 rT ^h fg - hL ]w 

dt = 3lT| 

dT t 
« 


(Mj-Wt) (4) 


dU 

3T 


abbreviating equation (4) into 
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au 

dT 


= lfl(T)l 1 V D + [f 2 (T)] 2 (M i -Wt) 


and substituting, gives an expression to be solved for the dewar saturation temperature as a function of 
time with given initial condition and constant thruster demand, W. 


• p Py n 

q-hfgW 1+ , T 

dT * [ 

dt [f,(T)l i V D + (f 2 (T)l 2 (M i -*t) 


Since the variables could not be separated, a numerical integration was required. The properties were 
defined by curve fits over the range of 1.4 to !.8°K. Methods and constants follow: 

q = Const = input - heat leak, J/sec (W) 

W * Const = input - demand flow rate, g/sec 
Tj * Const * Input - initial dewar temperature, °K 
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ORIGINAL PAGE Ml 

M{ ■ Const » input - initial helium mass, g ^ QUALITY 

Vq ■ Const ■ input - dewar volume, L 

Pj ■ Const * 145.1 g/L - helium liquid density, g/L 

Helium vapor is assumed to obey the perfect gas law 


p v * 


P 

RT 


g/L 


(P, tort; TK) ( 6 ) 


R* 15.5938 


torr L 

rr 


tew r _p_ i bp] g 

dT rx 2 * RT dT pK 


(P, torr; T, *K; 


3P tonr 
3T ,t K } 


(7) 


p . e < c o + c lO^) + c 2 t2 > 


torr (TfK) 


( 8 ) 


C 0 - 7.5959 
Cj - -9.9416 °K 
C 2 - 0.13913 V 2 

- I- qd/T 2 ) + 2 C 2 T] « (C ° + ° l(im + C2 l2) ton/°K OfK) (9) 

dT 


Using the Clausius Clapeyron relationship to get h fg 


d In P 
dT 



-_ + 2 c 2T 


h^ * (- C| R + 2 C 2 RT 3 ) Convi J/g 
3 hr- 

r ■ (6 C 2 RT 2 ) Convi J/g # K 
3T 


0?K) (10) 

CI?K) (11) 
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h L 
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dh L * C p dT 



Vl Coitvj dP 


■% 


Since P^ * Const, dv^/dT a 0 


P@T 


and 


f v L dP is negligible 


f 

t B 


h L * h L bMe *] C P dT 


with 


C p * c o + C 1 T + c 2 


h L()aje - 0. @ T baje - 1.6°K , then 

h L * c o T + ~Y T 2 + y- T 3 + D J/g (T°K) 

c Q = 8.1071 J/g°K 
cj *= -13.2199 J/g °K 2 
c 2 = 5.7071 J/g °K 3 
D = -3.842 J/g 


— » c 0 + cj T + c 2 T 2 J/g °K (T,°K) 


Convj - 0.13331 


J/L 

torr 


( 12 ) 


(13) 
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APPENDIX B 


HELIUM DELIVERY PRESSURE LOSS DERIVATION 




The pressure losses are a combination of frictional losses and expansion/contraction type losses 
and the loss for each section can be expressed as the sum of time two, 


AP= ft 


a/2 


W 


i/2 


+ K 


W 


D A 2 p2g c A 2 p2g c 


Since the flow is laminar, the friction factor, f, may be expressed as 


64. _ 64 

N Re = W D/A n 


and the average density, p, can be related to a nominal density, p Q , at a nominal pressure, P Q , and the 

average pressure, P, by assuming it obeys the perfect gas law P = pRT, and the temperature is constant 
across the resistance element. 


P = P 


o 


P 

K 


Combining and rearranging 


PAP 


64m LP 0 
A 2«cD 2 p 0 


w + 


Kp o 

A 2 2 *c «o 


w 2 


Applying this equation to a line represented by four sections as an example 
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P| APj -K,, W + K 21 W 2 
p 2 ap 2 - k 12 w + k 22 w 2 

P 3 AP 3 - Kj 3 W + K23 W 2 
P 4 ap 4 - k 14 W + k 24 W 2 

Summing 

(P, APj + P 2 AP 2 + P 3 AP 3 + P 4 AP 4 ) * (Kj j + K,2 + K| 3 + Kj 4 ) W + (K 21 + K 22 

+ K 23 + k 24 > - K 1T W + W 2 

The average pressure at each section is the upstream pressure at that section less half the pressure loss 
for that section v 

?i - ob - m 

f 2 * (P D - AP| - APj/2) 

P 3 ■ CP D - APj - APj - APj/2) 

?4 * ( P D - iP l - 4Pj - APj - AP4/2) 
and substitution into the above equation gives 



W + K 2 T W 2 -(p D AP 1 -^!^p d 


AP 2 - AP| AP 2 - 


ap 2 2 


-ap 2 ap 3 - 


AP3 2 


) + ( P D 

^ + ^P D AP 4 - APj AP 4 - AP 3 AP 4 - 


AP 3 - APj Ap 3 
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Kjt w + Kjt W 2 « P D UP] + APj + APj + AP^) - i (AP] ♦ APj ♦ AP3 + AP4) (AP] 
+ AP 2 + APj + AP^) 

Recognizing that the bracketed terms are the overall lots, AP t , the expression Rim p Hf frf to 
k 1t W + Kjt w 2 - p d ap t - i ap t 2 

Solving the quadratic in W gives 

. “ k 1T + >/ k 1T 2 “ 4 k 2T ( 7 aP T 2 “ P D APT) 

W ** £ 

2 ^ 

or alternately for aPj 

aP t * p d - \^d 2 - 2 (Ki T w + k 2T W 2 ) 
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HELIUM DELIVERY PRESSURE LOSS 
RESISTANCE FACTORS DETERMINATION 


Resistance Factors for Porous Plug 

Flow characteristics (Fig. S) for the porous plug were obtained from Reference 3. The sharp 
change of slope of these moves was referred to a the “critical mass flow rate** point and because of it the curve 
required two resistance representations. The unit area of the tested plug, Ap, was the same as that to be used on 

the GPB dewar. Temperature range for these equations is 1.6 <T< 1.8 °K. 

The dividing line for “ciidtal” flow was estimated to be 

*crit - U.fc + 10 (T - 1.65)1 x 10' 3 g/sec <T, °K) 

The pressure-flow relationship for flow below critical, i.e., low resistance, was 

AP * 0.5769 W torr (W, g/sec) 

P AP » 0.5769 P W . 

IfP*P D 

P AP » 0.5769 P D W torr 2 (W, g/sec; P D , torr) 


with 


K n « 0.5769 P D . 

Above critical flow, i.e., high resistance, 


W 


, L90 , ... 3.65 

1.65 + [T - 1.651 + - ■ 

°* 2 4.50 x 10 -3 


AP x 10- 3 g/sec (T, # K; AP, torr) 


AP * 1.233 W - 1.233 [1.65 x 1(T 3 (1 + 9.5 a/1.65 - 1))J 
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r s «2.i«K 


'TfllQK 


T a -1.«K 


>T§ ■ 1JMC 


CJ Pm IP 


Figure 5. Porous plug flow versus pressure drop. 


P AP « 1.233 P W - 1.233 P [1.65 x 1(T 3 (1 + 9.5 (T/1.65 - 1))1 


PAP - 1.233 P D W - 1.233 P D { 1.65 x 10‘ 3 (1 + 9.5 (T/1.65 - 1))1 torr 2 

-K n W+C (Pp.torrs T,°K) 


K| | - 1.233 P D 


C - -1.233 P D 1 1.65 x 1(T 3 (1 + 9.5 a/1.65 - 1))) . 
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Figure 6. Coding rings schematic. 
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OF POOR QUALITY 

The extra term C is not common to the other resistance dements since they do not have tike 
transition point in the W, AP curve. 

Equation for Helium Delivery System Pressure Loss Simulation for Computer Program 
~ g/sec 

AP and Pjj ~ torr 

AP = P D - J P D 2 - 2 (K 1t W + K 2T W 2 + C) 

W above W cr j t 

C = - 1.233 P D [1.65 x 10’ 3 (1 + 9.5 (Tjj/1.65 - 1))] 

K 1T = 121.97+ 1.233 P D 
K 2T =605.6 
W below ^ crit 
C = 0.0 

K 1T = 121.97 + 0.577 P d 
K 2 't a 605.6 

where 

W cr it — (1.5 + 10 (T - 1.65)) x 10 3 g/sec t 
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Figure 7. Distribution lines schematic. 
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p n 


iP 2 

ap 3 

^4 

* P S 

AP 6 

* P 7 


AC, 

iP 10 

^11 

iP K 

D 

R. 

R . 

IU 

R* 

IU 

R* 

IU 


IU 

Re 

IU 

Re 


* m * w 2 



ELE- 

MENT 

NO. 

ELEMENT 

L 

IN 

°h 

IN 

rf 

K c 

CONTR- 

ACTION 

*E 

EXPAN- 

SION 

*1 

J-? 

TORh- 

SEC/g 

*2> 

SttW 

TOW 2 

“p 

1 

POROUS PLUG BELOM CRIT. 

- 

- 

- 

- 

- 

0.5769Pq 


0. 

POROUS PLU6 ABOVE CRIT. 

- 

- 

- 

- 

- 

1.233 P 0 

- 

23 

■li 

2 

POROUS PLUG CONNECTOR 

6 

.335 

.0881 

- 

.68 

0.862 

60.8 


*ri 

3 

COOLING RING «1 

27.3 

.500 

2 * 

.2500 

- 

- 

6.861 

- 

*12 

4 

RING CONNECTOR 

8 

.460 

.1662 

.29 

.63 

7.146 

139.9 


5 

COOLING RING #2 

27.5 

.625 

2 * 

.3906 

- 

- 

16.881 

- 

*L3 

6 

RING CONNECTOR 

3.5 

.585 

.2688 

.40 

.66 

7.127 

185.0 

*R3 

7 

COOLING RING #3 

27.9 

.875 

2 * 

.7656 

- 

- 

11.633 

- 

R L4 

8 

RING CONNECTOR 

2.6 

.710 

.3959 

.55 

.65 

6.368 

174.4 

*R* 

9 

COOLING RING #4 

28.1 

1.000 

2 * 

1.000 

- 

- 

10.722 

- 

*ls 

10 

SUPPLY LINE 

10 

.96 

.724 

.28 

.47 

11.437 

42.6 

95 

1.71 

2.297 

- 

.29 

19.143 

2.3 

*SR 

11 

SUPPORT RING 

132 

1.25 

2*2* 

1.227 

- 

.33 

23.29 

.6 

“is 

12 

SPAR 

135 

3.9 

if/2 

11.95 



- 

.503 

- 


2 

BELOM CRITICAL FLOM 

121.97 

t0.5769(fc 

605.6 

° ! 

ABOVE CRITICAL FLOM 

121 .§7 
ft.233P n 

605.6 
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ORIGINAL p AGE EJ 
OF POOR QUALITY 


COMPUTER SIMULATION 


The results of the derivations and constants from Appendices A, B, and C were coded into the computer 
program shown in this section. 


Input 


NT 

Total number of time increments 

— 

NP 

Number of time increments between printouts 

— 

DTM 

Time increment for integration step 

Sec 

WI 

Initial helium mass in dewar (liquid + vapor) 

gram 

WDT 

Demand flow leaving dewar 

gram/sec 

Q 

Heat leak into dewar 

watts 

VD 

Dewar volume 

liters 

TK 

Upper temperature limit of helium in dewar 

°K 

TL 

Lower temperature limit of helium in dewar 

°K 

TI 

Initial temperature of helium in dewar 

°K 


Output 


Line 1 Input 


DTM 

WI 

WCT 

Q 

VD 

TH 

TL 

TI 

Time 

Initial 

Demand 

Heat 

Dewar 

Hi Temp 

Lo Temp 

Initial 

Increment 

Mass 

Flow 

Leak 

Vol 

Limit 


Temp 

sec 

gram 

G/sec 

Watts 

Liters 

°K 

°K 



Line 2 Initial Property Values 


PD 

RG 

PP 

PRV 

DH 

PDH 

H 

PHL 

Sat 

Press. 

P 

Sat Vap. 

Density 

Py 

3P/3T 


Latent 

Heat 

h fg 

3h fg /dT 

Liq. 

Enthalpy 

h L 

3h L /3T 

torr 

G/L 

torr/°K 


J/G 

J/G °K 

J/G 

J/G °K 


Line 3 Initial Pressure Loss 


PELP 

Delivery system pressure loss 

at initial conditions 

torr 
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Line 4 Initial Thruster Pressure 


PV 

Pressure delivered to thruster 

at initial conditions 

torr 


ORMJNAl PMC H 
Of POOR QUALITY 


Line 5 Current Property Values 


As line 2 without PD 

Property values at 

(current time - 1 time increment) 


Line 6 Current Dewar Conditions 


TM 

PD 

T 

DTDTM 

W 

Elapsed 

Sat 

Dewar 

Temp. Rate 

Helium Mass 

Time 

Press. 

Temp. 

of Change 

(Liq + Vap) 

t 

P 

T 

dT/dt 


sec 

torr 

°K 

°K/sec 

gram 


Successive Line Pairs 


Lines 5 & 6 Repeated 


Next to Last 


Same as Line 3 
for Final Conditions 


Last 


Same as Line 4 
for Final Conditions 
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C DEFINE PROPERTY FUNCTIONS OF POOR QUALITY 

PCT) - EXPCCO ♦ Cl/T ♦ C2*T**2) V 

RVCT) - P(T)/.(R*T> 

PPPTCT) - C2.*C2*T -.Cl/T**2)«P(T) 

PRVPT(T) ■ C-P(T)/T ♦ PPPTCT) >/(R*T) 

DHFO(T) C2.*C2*R*1**3 - C1*R)*CVI 
PDHPTCT) « CVt*6.*C2*R*T**2 

HUT) - D ♦ SCO*T ♦ (SCI/2. )*T**2 ♦ (SC2/3.)*T**3 
PHLPT(T) * SCO ♦ SCI*T ♦ SC2*T**2 
C 

READ(7 t 12) NT 
REAOC.7,12) HP 
READC J t 1 1 ) DTK 
READ.C 7, 1 1 ) Ml 
READC 7, II ) NOT 
READC 7 v 1 1 ) 0 
READ (.7* 1 1 ) VO 
REAOC7, II ) TH 
READC 7 • 1 1 ) TL 
READ (7, It ) T1 

NRITE(2« 10) DTM,HI.NDT,Q.VD,TH,TL,TI 
C SET UP CONSTANTS 
CV1 -.0.13331 
R - 15.594 
CO -.7.5959 
Cl - -9.9416 
C2 -.0.13913 
SCO > 8.1071 
SCI - -J3.2I99 
SC2. - 5.7071 
0 - -3.842 

C 

N-NT/NP 

C 

RL * 145.1 
C - INITIALIZE 
TM - 0.0 
T - TI 

C PRINT INITIAL PROPERTY VALUES 
PD ■- P(T) 

RC - RV(T) 

PP - PPPTCT) 

PRV - PRVPTCT) 

OH - OHFGCT) 

PDH - PDHPTCT) 

H * HLCT) 

PHL - PHLPTCT) 

NRITEC2,! 0) PD.RG, PP.PRV ,DH*PDH.H,PHL 
C CALC. PRESS. LOSS TO THRUSTER PLENUM 
PD - PCT) 

CALL PLOSSCNDT.PD.T) 
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C INTEGRATE TEMP* MATE. OF CHANGE FUNCTION ® 

DO 2 I- l f N OF POOR QUALITY 

DO I J-I.MP 
DH - DHFG(T) 

PRV -PRVPTCT) 

PDH ■ PDHPKT) 

RG - RV(T) 

PHL - PHLPT(T) 

RR - RG/CRL - RG) 

N -NI - WDT*TM 
IFtRG*VD.GE.N> GO TO 4 

Bl - VD*(DH*PRV*< I • * RR**2) ♦ RG*PDH*( 1 • ♦ RR) - CVI*PPPT(T)) 
82 * PHL - RR*PDH - DH*PRV*(I. ♦. RR)/(RL - RG) 

DTDTM * CO - DH*( I • ♦ RR)*NDT)/(BI ♦ 82*N) 

T - T * DTDTM*DTM 
PD ■ P<T) 

TM - TM ♦ DTN 
IF (T -TH> 3,3,4 

3 IF. <TL - T) 1,1,4 

1 CONTINUE 

NRlTE/2,.10) RG,PP,PRV,DH,PDH,H,PHL 

2 WRITE (2, 10) TM , PD, T, DTDTM, W 

4 CALL PLOSSI WDT,PD,T) 

10 F0RMAT(JX,9( I X, El 2.5)) 

11 FORM AT CFI 0.5) 

12 FORMAT! 16) 

STOP 

END 

C 

C 

SUBROUTINE PLOSS(WOT,PD,T) 

WDBO • WDT 

WDCRT* <1.5 ♦ IO.*(T — 1 .65))*! ,0E-3 
CK2T • 605.6 
IF ( WDBO WDCRT) 1,2,2 

C FLOW RATE BELOW CRITICAL 

1 C * 0.0 

CM IT ■ J21.97 ♦ <5.77E.-l)*PD 
GO TO 3 

C FLOW RATE ABOVE CRITICAL 

2 C- -I .233*PD*(1.65E-3*<». ♦ 9.5*XT/I.65 -.1.))) 

. CM IT." 121.97 4 I .2 33* PD 

3 DELP - PD - SQRT<PD**2 - 2.*(CMIT*WDB0 4. CK2T*WDBO**2 ♦ C)) 

PV * PD - DELP 

WRITE<2,43 DELP.Py 

4 FORMAT ( IM,E 1 2.5) 

. RETURN 

END 
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